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G6PDH is the rate-limiting enzyme of the pentose phosphate pathway and one of the principal source of NADPH, 
a major cellular reductant. Importantly, in ruminant's metabolism the aforementioned NADPH provided, is utilized for 
de novo fatty acid synthesis. Previous work of cloning the ovine {Ovis aries) og6pdh gene has revealed the presence of two 
cDNA transcripts (og6pda and og6pdb), og6pdb being a product of alternative splicing not similar to any other previously 
reported. 1 In the current study the effect of energy balance in the ovine G6PDH protein expression was investigated, 
shedding light on the biochemical features and potential physiological role of the 0G6PDB isoform. Changes in energy 
balance leads to protein expression changes in both transcripts, to the opposite direction and not in a proportional 
way. Negative energy balance was not in favor of the presence of any particular isoform, while both protein expression 
levels were not significantly different (P > 0.05). In contrast, at the transition point from negative to positive and on the 
positive energy balance, there is a significant increase of 0G6PDA compared with 0G6PDB protein expression (P < 0.001). 
Both 0G6PDH protein isoforms changed significantly toward the positive energy balance. 0G6PDA is escalating, while 
0G6PDB is falling, under the same stimulus (positive energy balance alteration). This change is also positively associated 
with increasing levels in enzyme activity, 4 weeks post-weaning in ewes' adipose tissue. Furthermore, regression analysis 
clearly demonstrated the linear correlation of both proteins in response to the WPW, while energy balance, enzyme 
activity, and 0G6PDA relative protein expression follow the same escalating trend; in contrast, 0G6PDB relative protein 
expression falls in time, similar to both transcripts accumulation pattern, as reported previously. 2 



Introduction 

Studying the biology of adipose tissue in farm animals and 
acquiring the knowledge of animal production and health, we 
may comprehend better the significant resemblances in lipid 
physiology with emphasis in lipogenesis. There are species-based 
biases that occur, not only in adipose tissue features but also in 
the metabolism, as a whole. In all species studied so far, g6pdh 
is considered as an X-linked gene displaying a "housekeeping" 
profile, while several variants exist. 1,2 Having mentioned these 
biases in the metabolic pathways, lipogenesis unfolds the hidden 
selective pressure that nature reveals by the adaptive mechanisms 
that led ruminants to divert from other herbivorous mammals in 
order to survive, share mutually natural resources and interact 
with the other species, as well. 

Focusing on the aforementioned metabolic pathway and spe- 
cifically the biosynthesis of fatty acids starting from acetyl-CoA, 
there are certainly considerable amounts of reducing agents, 
like NADPH, which are required for the reductive pathway and 
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leading to fatty acid synthesis. Major important enzymes linked 
to NADPH production are: cytosolic isocitrate dehydrogenase 
(IDH), EC 1.1.1.42, cytosolic NADP malic enzyme, ECU. 1.40, 
enzyme of pentose phosphate shunt, glucose- 6 -phosphate dehy- 
drogenase (G6PDH), EC 1.1.1.49, and 6-phosphogluconate 
dehydrogenase, (6PGDH), EC1. 1.1.44. The last two enzymes are 
the main NADPH suppliers in ruminants, providing 50-80% of 
the total required NADPH for the fatty acid synthesis. 3 

G6PDH is a cytosolic enzyme encoded by a housekeeping 
X-linked gene whose main function is to produce NADPH, a 
key electron donor in the defense against oxidizing agents and in 
reductive biosynthetic reactions. It is conventionally referred to, 
as the rate limiting enzyme of the pentose phosphate pathway in 
all cells, 4 and it is therefore regarded as important in the biosyn- 
thesis of the sugar moiety of nucleic acids. In addition, G6PDH 
apart from its role as NADPH "generator", the control of key 
enzymes modulating tissue response, such as llbeta-hydroxys- 
teroid dehydrogenase-1 (11(3-HSD-1) that is highly expressed 
in adipose tissue, liver, and brain, is an additional role. 5 In the 
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aforementioned tissue, G6PDH acts primarily as an oxo-reduc- 
tase to generate active glucocorticoid placed in the endoplasmic 
reticulum (ER). 5 

As far as the lipogenic activity of G6PDH is concerned, it has 
been demonstrated that hepatic G6PDH is regulated by nutri- 
tional signals, including high-carbohydrate diet, polyunsaturated 
fatty acids, and hormonal signals such as insulin, glucagon, thy- 
roid, and glucocorticoids but is also participating in the reduc- 
tive biosynthesis of fatty acids and cholesterol. 6 Previous research 
reports that enzymatic activity and expression levels of G6PDH 
were significantly elevated in, with pleiotropic role, white adipose 
tissue of obese models, and diet-induced obesity mice. 7 It is well 
known in obesity, that adipose tissue exhibits increased oxidative 
stress, which is a major contributor to metabolic disorders such as 
insulin resistance, playing also a central role in nutrient sensing 
and energy homeostasis. 

Concerning ovine g6pdh, several studies have been conducted 
on enzymatic level showing the reaction of enzyme activity to 
different stimuli. 8 " 10 At the molecular level, g6pdh has been well 
studied in many eutherian species. The g6pdb gene, is remark- 
able for its genetic diversity and there are many variants of 
g6pdh, mostly produced from missense mutations, that have 
been described with wide range levels of enzyme activity, while 
in humans are associated with clinical symptoms. 11 Complete 
characterized sequences have been already published for human, 
mouse and rat genes, 12 " 15 while many studies have shown the reg- 
ulation of gene expression in respect to various stimuli. 2 G6PDH 
is an ubiquitously expressed molecule and glucose-6-phosphate 
dehydrogenase enzyme activity is in accordance with the energy 
uptake. 2 

Previous work in our laboratory, has reported the ovine g6pdh 
gene cloning and characterization, 16 " 18 elucidating the presence 
of two differentially expressed ovine g6pdh transcripts (og6pda 
and og6pdb). The og6pdb gene product is the result of a splicing 
event leading to a transcript with 27 extra nucleotides, in total. 
This extra segment causes a frameshift in the polypeptide chain 
resulting in changes around the area of the substrate binding site 
as an enlarged binding "pocket". It has been hypothesized that 
the 0G6PDB transcript could probably result in a low G6PDH 
activity, 16 as analogous events have been previously reported in 
sheep. 19,20 Thus, the presence of a novel alternative transcript of 
ovine g6pdh gene, mirrors a potentially novel role in lipid metab- 
olism. Another published work 21 from our lab has studied bio- 
chemical parameters for metabolic enzymes involved in lipogene- 
sis in ewes from the same breed and at a similar lactation stage, as 
in this work. That report has shown the effects of energy balance 
on a few biochemical parameters like insulin, free fatty acids, and 
(3-hydroxybutyric acid concentrations. 21 In detail, lower values 
for insulin, higher levels of growth hormone (GH) and increased 
rates for the latter two acids were observed at the negative energy 
balance in the period of the first 2 weeks of lactation, while this 
period lasted only 5 days in Karagouniko breed. 21 These results 
share similarities with this work, concerning the negative energy 
balance "dominant" protein isoform, 0G6PDB. As in our work, 
we had intense interest to investigate the mechanism of og6pdh 
gene regulation upon energy balance changes. We conducted our 
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Figure 1. Representative results of immunoblot analysis of protein 
extracts of tail-biopsies of animals fed on (1) negative; (2) transit to posi- 
tive; and (3) positive energy balance. Twenty-five micrograms of protein 
were loaded on a 8-16% gradient polyacrylamide gel and dry-blotted 
on PVDF membrane and further incubated with primary G6PDH-specific 
and HRP-secondary antibody (see Materials and Methods). Samples 
numbers correspond to the following samples: (1) first WPW; (2) second 
WPW; and (3) fourth WPW. Protein expression of 0G6PDA and 0G6PDB 
was altered through the course of time and the levels of p-actin protein 
are shown below in the immunoblot. 0G6PDH protein isoform expres- 
sion was normalized to fS-actin protein. 



research at the translational level, to identify the putative protein 
0G6PDH isoforms. 

In order to explore the above role, our research deal with the 
effect of dietary energy level in ovine G6PDH protein isoform 
expression, to elucidate the presence of the two ovine G6PDH 
protein isoforms, the pattern of expression of both protein iso- 
forms as influenced by the energy balance in the course of time 
post weaning, and the potential physiological role of og6pda/b 
gene products. Evidence from other species for G6PDH variants 
in electrophoretically separated protein extracts 22,23 was the aspi- 
ration in our research, while the pertinent investigation line was 
lacking in ruminants since the discovery of the second og6pdb 
transcript. 16 Our findings not only shed light and confirm simi- 
larities and differences in both, og6pdh transcript accumula- 
tion 24 and 0G6PDH isoform protein expression, in respect to the 
energy balance changes, but also emphasize the significance of 
the alternative splicing g6pdb product in the regulation of g6pdh 
gene expression and 0G6PDH enzyme activity. 

Results 

Results from our immunoblot analysis shown in Figure 1, 
were further analyzed by quantitative densitometry and gener- 
ated the mean values in Figure 2A and B. In our Figure 2A graph 
we summarize the data from 0G6PDA and 0G6PDB protein 
expression variations as a function of the energy balance in ewes' 
adipose tissue. As reported elsewhere by the authors, 24 the average 
milk production as a yield was not statistically different among 
the week groups post-weaning. Nevertheless, the energy balance 
changed at the end of first, second, and fourth WPW, i.e., -0.50 
± 0.24 MJ, -0.19 ± 0.14 MJ, and 4.48 ± 0.26 MJ respectively, 24 
with significant increase toward the fourth WPW. 
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Figure 2. Average relative expression of the 0G6PDH protein isoforms in 
the ovine adipose tissue. (A) Average relative expression of the 0G6PDA 
and 0G6PDB protein. Relative amounts of protein were determined in 
each sample with (3-actin protein as the housekeeping gene product 
and endogenous control. Each sample was analyzed from each biopsy 
consisted of 7 animals. (B) Average relative expression of 0G6PDA and 
0G6PDB proteins in the ovine adipose tissue. Data were merged in 2 
groups according to their negative or positive energy balance. 0G6PDA 
and 0G6PDB protein expression vs. negative or positive energy bal- 
ance. Relative amounts of protein were determined in each sample and 
(3-actin protein as the endogenous control. Values are in arbitrary units 
(mean ± SEM, n = 6,P> 0.05, ns, not significant). 



Our densitometry data in Figures 2A and 3A show the rela- 
tive protein expression of 0G6PDA and 0G6PDB isoform, where 
energy balance switch from the negative (first and second week) 
to the positive (fourth week) giving rise to a concomitant trend 
in protein expression of G6PDA. Specifically, 0G6PDA protein 
isoform rises from 1.052 ± 0.071 (first week) and 1.265 ± 0.209 
(second week) to 1.596 ± 0.1982 (fourth week) with significant 
changes between first and fourth week (P < 0.05). On the con- 
trary, 0G6PDB relative protein expression shows a more signifi- 
cant decline which is observed with more significant changes 
(P< 0.01) between the first and second WPW from 1.011 ±0.096 
and 0.548 ± 0.099 respectively, to the fourth week (0.528 ± 0.058), 
compared with the 0G6PDA relative protein expression levels. 



Importantly, the calculated 0G6PDA/0G6PDB ratio from 
Figure 2A in protein expression between these two protein iso- 
forms shows an increase in the course of time (P< 0.01, P < 0.001) 
from the first week (1.077 ± 0.0610) and second week (2.515 ± 
0.467) to the fourth week (2.977 ± 0.218) (data not shown), which 
is consistent with the similar pattern of the g6pda/g6pdb tran- 
script accumulation ratio, reported elsewhere. 24 Hence, there is a 
significant increase in the 0G6PDA compared with the 0G6PDB 
protein levels, and the attributed 0G6PDA/0G6PDB ratio (data 
not shown) in the protein expression level increase is due to the 
smaller scale changes in G6PDB protein expression, compared 
with 0G6PDA (Figs. 2A and 3A). Specifically, there are signifi- 
cant changes in both protein isoform expression focused, on one 
hand in 0G6PDA protein isoform expression changes between the 
distal time points i.e., first and fourth WPW, while on the other 
hand, statistical analysis concerning 0G6PDB protein isoform 
expression shows significant changes between the first and second 
as well as between first and fourth WPW. Interestingly, signifi- 
cant changes in both protein isoform expression were not observed 
between the second and fourth WPW (Figs. 2A and 3A). 

By organizing our data according to the energy balance (abso- 
lute negative or positive values) in 2 groups, ewes in negative 
energy balance and ewes in positive energy balance, and follow- 
ing similar statistical analysis the graph in Figure 2B, was plot- 
ted, illustrating relative protein expression values versus the group 
of either negative or positive energy balance. This "intervention" 
was useful to understand the "net effect" of energy balance in 
WPW. It is clearly shown, that negative energy balance 0G6PDA 
protein expression had an almost 2-fold increase (P < 0.01) com- 
pared with 0G6PDB, i.e., (1.294 ± 0.143) compared with (0.678 
± 0.088) respectively. Equally, in the group of the positive energy 
balance 0G6PDA protein expression had at least 2-fold increase 
(P < 0.01) compared with 0G6PDB, i.e., (1.332 ± 0.129) com- 
pared with (0.678 ± 0.088) respectively. Nevertheless, the energy 
balance change did not affect significantly, either 0G6PDA or 
0G6PDB relative protein expression (P > 0.05) (Fig. 2B). Fur- 
thermore, the total contribution from both isoforms in 0G6PD 
protein expression from the negative energy balance data group, 
to the positive, from first WPW to fourth WPW results in a 
fall at the protein level, nevertheless these changes between the 
negative and positive energy balance groups were not statistically 
significant for both protein expression (0G6PDA and 0G6PDB), 
and only affected expression changes between the isoforms, as 
shown from the graph (Fig. 2B). 

To further analyze the energy balance effect on protein 
expression levels, regression analysis of the data was performed as 
shown in Figure 4A and B. Both 0G6PDA and 0G6PDB relative 
protein expression linear regression analysis, demonstrate how 
far protein expression and energy balance values deviate from a 
putative linear relationship. The R squared values were estimated, 
7? oG6PDA 2 = 0.25 and ^ oG6PDB 2 = 0.35 and post-hoc analysis test and 
best-fit values estimates for the y intercepts at x = 0.0, were a oG6 
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< 0.05, P oG6mB < 0.01). Interpreting these data, is clearly 
demonstrated the linear correlation of both proteins in 
response to the WPW and/or energy balance changes. 
In Figure 4A and B, is also shown from the graphs that 
the cross-section of both lines characterizes no change 
in protein isoform expression (equilibrium state), dur- 
ing the biopsy week and/or energy balance (MJ), at the 
protein expression level. Interestingly, from the graphs 
(Fig. 4A and B), this cross-point is located at the rela- 
tive protein expression level, ca. 1.0. Under these condi- 
tions and at approximate protein expression level, ca. 
0.96, this equilibrium occurred before the end of first 
biopsy week (Fig. 4B). The same cross-point also coin- 
cides with negative energy balance, where both lines 
meet at ca. -6.5 MJ (Fig. 4A). 

As far as the 0G6PDH enzyme activity is concerned, 
Figure 3A shows statistically significant values increase 
toward the fourth WPW, where this activity reaches 
maximum values on the fourth WPW, as shown by the 
enzyme activity, nanokatal (nkat) at the right jy-axis, in 
Figure 3A. Figure 3A also shows a comparison with 
both protein isoform expression levels to demonstrate 
the concomitant trend of 0G6PDA and the contrary 
trend of 0G6PDB relative protein expression, simulta- 
neously with the escalating G6PDH enzyme activity in 
the samples. 

In order to evaluate our data further and rule out 
issues, such as differences in increased enzyme content 
vs. activity, we plotted on a graph relative expression 
as a function of enzyme activity (Fig. 3B). As shown 
in Figure 3B, enzyme activity starts from a minimum 
(ca. 6.93 nkat, first WPW), at equilibrium point where 
both isoforms are equally expressed (cross-point M in 
Figure 3B), reaching a maximum in the enzyme activ- 
ity, which has been detected at, ca. 62.06 nkat, fourth 
WPW) for both protein isoforms. At this point, where 
the maximum of values in 0G6PDA relative protein 
expression is 1.596 ± 0.198, compared with the lowest 
values in 0G6PDB relative protein expression which is 
0.528 ± 0.058, in the fourth WPW, the relative pro- 
tein expression per enzyme activity is 0.026 ± 0.023 
and 0.008 ± 0.003 (Fig. 3B). A positive association is 
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Figure 3. (A) Left-side graphs demonstrate the average relative expression of 
0G6PDA and 0G6PDB proteins in the ovine adipose tissue, compared with the 
right-side graph of 0G6PDH enzyme-specific activity (nkat) vs. WPW. One way 
ANOVA vs first WPW was performed to evaluate significant changes followed by 
the Bonferroni multiple comparison test. Values are in arbitrary units (mean ± SEM, 
*P < 0.05, **P < 0.01, ns = not significant). (B) Average relative expression of the 
0G6PDH protein isoforms in the ovine adipose tissue vs 0G6PDH enzyme activity. 
x-y graph was plotted from mean ± SEM in every point. Points I, II, and IV indicate 
first, second WPW, and fourth WPW, respectively. Points II and IV show statistically 
significant differences between the two isoforms. M symbol is the cross-section of 
the two curves, where 0G6PDA and 0G6PDB has equal relative protein expression, 
and enzyme activity starts at 7.256 ± 3.040 (nkat). The values from the curve are rel- 
ative expression (arbitrary units/nkat) (mean ± SEM, n = 6, **P < 0.01, ***P < 0.001). 
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(Fig. 3B), for the 0G6PDH enzyme activity changes (nkat) 
(x-axis), with enzyme activity starting from 7.256 ± 3.040, 
17.82 ± 5.027, with maximum at 62.06 ± 14.985 nkat (P < 
0.01), at time points I (first WPW), II (second WPW), and IV 
(fourth WPW), respectively (Fig. 3B), and 0G6PDA relative 
protein expression starting from 1.052 ± 0.071, 1.265 ± 0.210 
with maximum at 1.596 ± 0.198 arbitrary units, at the same 
time points (Fig. 3A and B). This is in contrast with 0G6PDB 
relative protein expression (opposite trend) starting from 1.011 
± 0.096, 0.549 ± 0.099 with minimum at 0.528 ± 0.058, from 
the negative to the positive energy balance (first, second, and 
the fourth WPW, respectively) (Fig. 3A and B). In both cases, 
the linear trend for G6PDA expression and enzyme correlation 
is met with good correlation coefficient (r 2 = 0.954) in the range 



of enzyme activity values, whereas correlation coefficient for 
0G6PDB expression and enzyme activity is lower (r 2 = 0.459). 
In detail, 0G6PDA protein isoform expression continues to rise 
in a linear fashion with the enzyme activity (r 2 = 0.954), while 
0G6PDB protein isoform expression continues to fall, reaching 
an almost plateau (value at ca. 0.008 ± 0.003 expression/nkat) 
(Fig. 3B). Interestingly, further statistical analysis has shown, 
that the former curve fits well to a classic linear curve while the 
latter fits better to one-phase decay curve compared to the linear 
form estimated earlier (i.e., goodness of fit, r ,. 2 = 0.57 > 

v o ' nonlinear 

r lm J = 0.459) (Fig. 3B). 

So conclusively, we observed changes with the same trend, as 
has been shown previously for g6pda relative transcript accumu- 
lation 24 and 0G6PDH protein expression vs. enzyme activity, as 
well. Figure 3B, is also emphasizing the fact described earlier in 
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Figure 4. (A) Protein expression ratio (relative to (3-actin) of 0G6PDH 
isoforms (0G6PDA and 0G6PDB) in ovine adipose tissue as a function 
of WPW. Linear regression analysis was performed to test the linear 
relationship between protein expression of protein forms of 0G6PDH 
and energy balance. A positive 0G6PDA and negative 0G6PDB protein 
expression correlation was found in relation to the energy balance 
changes. The slope differences between the best-fit curves were sig- 
nificant (***p< 0.001). The insert depicts the equation formulae for both 
best-fit lines. (B) Protein expression ratio (relative to (3-actin) of 0G6PDH 
isoforms (0G6PDA and 0G6PDB) in ovine adipose tissue as a function of 
energy balance (MJ). Linear regression analysis was performed to test 
the linear relationship between protein expression of protein forms of 
0G6PDH and energy balance. A positive 0G6PDA and negative 0G6PDB 
protein expression correlation was found in relation to the energy bal- 
ance changes. The slope differences between the best-fit curves were 
significant (***p < 0.001). The insert depicts the equation formulae for 
both best-fit lines. 



Figures 2A and 3A that significant changes in relative protein 
expression were observed at the time points (WPW) where the 
detected enzyme activity is in the second or fourth WPW. The 
significant changes between the two 0G6PDA and 0G6PDB 
isoform protein expression, were observed between the second 
and fourth WPW, and not at point I (first WPW) where enzyme 
activity at its lowest value (7.256 ± 3.040) (Fig. 3B). Interestingly, 
the distance of these expression curves after the second WPW, 
tend to a parallel direction reaching their maximum difference in 



relative protein isoform expression with the same enzyme activ- 
ity (ca. 62 nkat), i.e. maximum for 0G6PDA and minimum for 
0G6PDB (1.596 ± 0.198 vs 0.528 ± 0.0585 respectively). 

Discussion 

Previous studies have reported the cloning and molecular 
characterization of the ovine g6pdh gene as well as its promoter 
region. 16 " 18 These studies have analyzed and confirmed the pres- 
ence of two og6pdh transcripts at the cDNA level. These alter- 
native spliced og6pdh gene products, i.e., the expressed protein 
isoforms as shown in our data, accumulate evidence on the dis- 
parate physiological role of these transcripts. By studying the 
influence of the energy balance alterations on the og6pdh tran- 
script accumulation 24 and protein levels, these patterns of expres- 
sion could shed light on the translational control of the 0G6PD 
protein product(s) and further elucidate the potential physi- 
ological role of the recently identified second (g6pdb) alternative 
transcript. 16,24 

As, shown from our protein analysis and immunodetection 
results, on both transcript products i.e., 0G6PDA and 0G6PDB, 
there is a concomitant rise in the accumulation of these pro- 
teins on the first, and a falling pattern on the second protein 
respectively, following the transition from negative to positive 
energy balance, in the adipose tissue of the ewes. This change is 
explained from the fact, that animals fed with high carbohydrate 
diet are led to higher energy balance, i.e., high energy excess and 
positive energy balance, higher enzyme activities in the lipogen- 
esis, such as fatty acid synthase, acetyl CoA carboxylase, malic 
enzyme, and 0G6PDH expression levels, compared with the 
previous lower energy state. All the aforementioned factors are 
coordinated depending on the energy demands at a certain time, 
resulting to the enhancement or attenuation of the final g6pdh 
gene product and finally its protein synthesis. Previous work 
regarding the aforementioned changes at the g6pdh mRNA and 
protein level has been extensively reported. 2,25 

Published work has shown that G6PDH-deficient mice had 
increased oxidative stress, decreased NADPH levels, decreased 
glutathione reduced form (GSH) levels, as well as increased 
markers of lipid peroxidation. 26 In the above oxidative reactions 
glutathione contribution is rate-limiting depending on cysteine, 
which is the rate-limiting factor in cellular glutathione synthe- 
sis. 27 G6PDH is also essential for embryo development, as the 
complete knockout of g6pdh gene was embryonically lethal. 28,2 ' 

In our data, the observed values i.e., the rising trend in protein 
expression of 0G6PDA from the negative to the positive energy 
balance, is concomitant with the other estimated parameters. 
Hence, is directly related to 0G6PDA transcription level, 24 pro- 
tein expression and total enzyme activity (Figs. 2A and 3A and 
B). In addition to that, og6pdb alternative transcript accumula- 
tion, 24 as well as protein expression, show similar but opposite to 



da transcript trend (Figs. 2A and 3A and B); however the 
0G6PDA protein preponderance over 0G6PDB protein in the 
course of the WPW is apparent. 

Previous work in our lab has reported another study for the 
g6pda/b transcripts, 24 which in comparison with our data at the 
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protein level, suggest that og6pdh gene transcript accumulation 
is more abundant and stable, in terms of producing mRNA, 
than total 0G6PDH protein, while this is also more intense for 
0G6PDB isoform. 24 This difference with the previous report 
is also apparent from our Figures 2A and 3A, where the max 
x-fold levels of protein expressed were below 2x-fold for 0G6PDA 
and almost 0.5^-fold for 0G6PDB protein expression during the 
fourth WPW, compared with 3x-fold and 1.5x-fold expression 
in relative transcript levels for og6pdh gene respectively. 24 There- 
fore, estimating from the above observations, the relative pro- 
tein expression to the transcript accumulation ratio is around 
-66% for 0G6PDA and -33% for 0G6PDB, suggesting that 
the level of expressed proteins are marginally correlated with the 
level of mRNA transcripts. This is not unusual, as changes in 
gene expression level are frequently not reflected at the protein 
level. The current gene regulatory events are not known and 
require further investigation, nevertheless both transcripts/pro- 
teins share equal importance. This intensive transcription yield- 
ing lower 0G6PDH isoform protein levels could be investigated 
further with more promising scrutiny, by microarray technology 
with studies of functional regulation that could, in principle, be 
extended to detect not only the changes in the overall expres- 
sion of a gene, but also changes in its splicing pattern between 
different tissues. Interestingly, the results from the translational 
events in these two transcripts reveal a differential regulation, as 
responses and roles between them in adipose cells apparently dif- 
fer. In fact, alternative splicing events consist of a control mecha- 
nism in gene regulation, enzyme activity, and function, which 
generally adapt the functional enzyme production according 
to the metabolic conditions, in time. It is very likely, that this 
is applied in 0G6PDH protein expression and enzyme activity 
where the functional 0G6PDH protein is reflected on the ratio of 
these isoforms, as mentioned earlier. 

As shown from the graphs (Figs. 2A and 3A and B), the 
2 og6pdh gene transcripts and its products are functioning in 
two opposite directions to supply the cell with the required 
0G6PDH protein isoform levels and enzyme activity. Moreover, 
the ratio obtained by calculating the relative protein expression 
of the two isoforms, i.e., 0G6PDA/0G6PDB ratio (data not 
shown), show an upward trend from the negative to the positive 
energy balance, due to the greater protein level changes in favor 
of 0G6PDA. These data become more apparent as presented in 
Figure 3B. By this graph, we interested to show what the cor- 
relation is between protein expression for these two protein iso- 
forms, as detected by the immunoblots, and the enzyme activ- 
ity, in order to rule out differences in increased enzyme content 
vs. activity. As shown from Figure 3B, it is possible to have no 
change in enzyme content (plateau area in these curves), but 
measurable changes in activity, as the number of protein mol- 
ecules remain invariable, but enzyme activity reaches a maxi- 
mum, modulated possibly by negative feedback of the isoforms 
themselves (Fig. 3A and B). Importantly, 0G6PDB isoform 
expression is the lowest in terms of protein expression for the 
maximum measured enzyme activity in return. It seems that at 
the positive energy balance until the fourth WPW, less 0G6PDB 
protein is required in order 0G6PDH enzyme activity remains 



unaffected fulfilling its role in the adipose cell, simultaneously 
with 0G6PDA protein isoform (Fig. 3B), We cannot estimate 
activation/inhibitory sites on these 0G6PDH enzyme activi- 
ties and the possible mechanism at this stage, or the enzyme 
localization in terms of physiological significance. Although, we 
may assume that the changes in enzyme activity (possibly due 
to, either increase in the number of activation sites in 0G6PDA 
isoform or blend with molecules of 0G6PDB as this could result 
in lower 0G6PDH activity 16,19,20 ) could be also due to post- 
translational modifications, and the dominant protein isoform 
(0G6PDA) at the maximum level, is "rendered" by these modi- 
fications more active and remain possibly the major 0G6PDH 
activated enzyme at the positive energy balance. Interestingly, 
the estimated intercepts of both protein expression (£ oGSPDA 
and ^ oG6PDB ) show, that expressing 0G6PDA protein there is a 
ca. 1.30x-fold more direct response, as an effect of the energy 
balance changes, than 0G6PDB protein isoform (Fig. 4). This 
point will be further elucidated by experimenting on the indi- 
vidual enzyme affinities of the two isoforms for their substrates, 
or ideally by producing their recombinants and further scruti- 
nize their kinetic parameters. 

0G6PDH enzyme regulation and activity along with the 
energy balance, is functionally linked to the protein compo- 
nents and expression levels of the total 0G6PDH in terms of 
protein and enzyme activity. Although, we do not know the 
mechanism in this activity and the link with these protein com- 
ponents, nevertheless this is reminiscent of a balanced protein 
regulation from a more or less "stable" pool of the 0G6PDH 
protein from various deposited resources, reserved to be chan- 
neled to produce the enzyme and its activity on demand of 
energy requirement. 30 Thus, at low glucose, or energy, the flux 
of G6PDH for ATP generation to the cytosol, with cycling, 
is adjusted by the protein levels and the isoform in favor, in 
order to maintain a threshold redox state of NADPH. When 
higher levels of energy is provided through the diet, there is an 
increased flux through cytosolic 0G6PDH enzyme, while adi- 
pose tissue provides sufficient G6PDH for both ATP synthesis 
rates and cytosolic NADPH generation, which is allosterically 
stimulated by NAPD* to increased rates of ATP and NADPH 
synthesis. 30 As shown in Figure 4, the equilibrium point 
(cross-point), where both 0G6PDA and 0G6PDB were equally 
expressed in the first WPW (ca. third day post-weaning, the 
relative expression of the two isoforms appeared at the negative 
values of the curve and energy balance (tick point -6.69 MJ), 
which is interpreted as less energy balance input and consump- 
tion than what is required. Therefore, under these conditions 
and taking into account the almost linear relationship (Fig. 4), 
these results imply a starvation regime for the organism, where 
beyond the cross-point and toward the negative energy balance 
values, an active role is undertaken by 0G6PDB. We cannot 
predict at this energy stage, if 0G6PDB is a potential "regula- 
tor" of the other protein isoform (0G6PDA) or if both are regu- 
lated by the same control mechanism. Under this mechanism, a 
negative control feature emerges in 0G6PDH activity to protect 
homeostasis, restoring the necessary enzyme activity and keep 
an equilibrium between lipogenetic and lipolytic processes, as 
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Figure 5. Theoretical model of the role of 0G6PDH transcripts in energy metabolism. The picture depicts 
the influence of energy balance on ovine G6PDH mRNA and protein isoforms. Negative energy balance has 
inversely proportional effect on 0G6PDA mRNA accumulation and protein isoform. In contrast, positive 
energy balance has a positive effect on 0G6PDA mRNA and negative effect on 0G6PDB mRNA accumula- 
tion and protein isoform. (1) State of the cell in tissues producing and secreting both transcripts into the 
cytoplasm. The mRNA molecules produced after processing and alternative splicing, export to the cyto- 
sol and translation follows. Each isoform production, depends on the (-) or (+) energy balance. (2 and 3) 
Negative and positive energy balance favors either 0G6PDB or 0G6PDA mRNA and protein isoform. The 
molecules of 0G6PDH (° or ") as depicted in the circles (2 and 3), are for simplicity only (see Discussion). 
Thick lines and signs show the favored transcript/protein isoform production upon the energy balance 
state; ptm, posttranslational modifications; ± AG, ± energy balance. 



dictated by the energy balance. The energy balance dependence, 
also reveals a fine tuning feature with a more sophisticated and 
complex functional control of G6PDH enzyme (s). 

In the model that shaped our thoughts, we show graphically 
this interaction in detail, as depicted in Figure 5, for the role 
of 0G6PDH transcripts in energy metabolism. Low energy bal- 
ance diminishes og6pda transcript production in favor of og6pdb. 
The same trend is followed for the respective protein isoforms 
(0G6PDA and 0G6PDB). There is a contrary effect exerted by 
the positive energy balance on the transcript accumulation and 
protein isoforms, i.e., this favors og6pda gene product. Briefly, 
upon energy balance changes, more molecules are recruited from 
the favored protein isoform to supply the enzyme demands, at 
one time (see circles 2 and 3 in Fig. 5). Furthermore, the og6pdh 
gene products seem to have more than one roles to serve, briefly, 
on one hand in providing NADPH for cytosolic synthesis of fatty 
acids etc., and on the other hand to regulate key enzymes modu- 
lating responses such as llbeta-hydroxysteroid dehydrogenase 
1 (11|3-HSD1) in adipose tissue, localized in the endoplasmic 



reticulum (ER) besides 0G6PDH, 
as shown in Figure 5, revealing a 
multi-functional profile for this 
metabolic enzyme. 

In the disease front, diabetes 
research has shown that inhibi- 
tion of G6PDH can be used to 
inhibit insulin resistance in obese 
patients who developed insulin 
resistance. 31 They are several iso- 
forms of human G6PDH, but 
which variant is targeted/inhib- 
ited in the insulin resistance, or 
other existing chronic diseases 
including heart failure and pul- 
monary hypertension, has yet to 
be elucidated. 32 Under oxidative 
stress conditions (polyol path- 
way), the regeneration of reduced 
GSH is compromised. 33 G6PDH- 
deficient mice have increased 
oxidative stress and decreased 
GSH levels, as the reduced lev- 
els of NADPH are deficient to 
maintain the level of glutathione 
in the tissues, that protect them 
against oxidative damage. 34 Thus, 
G6PDH enzyme is regulated by 
nutritional and hormonal factors, 
at posttranslational and posttran- 
scriptional level as well as oxidant 
stress. 34 It would be interesting to 
know what is the outcome of the 
inhibitory effect, i.e., by transla- 
tional blocking or gene "knock- 
out" technology, on one or the other og6pdh gene product iso- 
forms. Moreover, the question to be addressed at this stage is 
what the control mechanism is that "senses" the energy balance 
changes prior to gene transcription/splicing/protein isoform 
"adaptation", adjusting the 0G6PDH enzyme activity. 

In conclusion, that transcript accumulation as reported 
before 24 is in agreement with the relative protein expression of 
both isoforms, nevertheless there is difference in the scale and 
intensity of changes between the transcript and the protein levels, 
where in the latter the level of signal is lower. However, all values 
reproducibly show similar trend in the relative 0G6PDH tran- 
script variants as detected by real-time PCR and relative protein 
isoform expression, by immunoblotting. These energy-depen- 
dent variables of og6pdh gene products, involved in lipogenesis/ 
lipolysis, are adaptable but concurrent with the energy balance 
changes, manifesting similarities at two levels of control — tran- 
scription and translation — as well as differences in scale of their 
response, in the course of time (weeks post weaning), in the adi- 
pose tissue of the ewes. 
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Materials and Methods 

General (animals, housing, and feeding) 

Seven randomly selected lactating ewes of the Chios breed 
were employed for the study. The animals were in the middle 
stage of lactation and housed at the experimental farm of the 
Agricultural University of Athens. The experimental design has 
been described in detail, in a previous report 24 ; briefly: ewes were 
fed individually twice per day, at about 7:00 and 16:00 h on an 
alfalfa hay (4.14 MJ/kg) and on a concentrated diet. The concen- 
trated diet consists of a basal (7.1 MJ/kg) and a lactation (7.3 Ml/ 
kg) ration. By the end of the second experimental week the 
amount of diet was set in such level that a negative energy balance 
was achieved according to ewe's requirements (0.4 kg alfalfa hay, 
0.5 kg basal ratio, 1 kg lactation ratio). The next two weeks, the 
amount of diet achieved a positive energy balance (0.4 kg alfalfa 
hay, 0.5 kg basal ratio, 1.5 kg lactation ratio). Water was offered 
free ad libitum under all experimental conditions. At the end of 
the first, second, and fourth experimental weeks, a tail biopsy 
was operated on all animal groups from the subcutaneous adi- 
pose tissue. Prior to biopsy ewes were fasted but had free water ad 
libitum. At the day of sampling ewes were anesthetized with the 
use of Tiletamine and Zolazepam (Zoletin 50) in 15-20 mg/kg 
BW. Samples of adipose tissues were snap frozen in liquid nitro- 
gen and stored at -80 °C for further protein extraction. All exper- 
imental procedures were approved by the bioethical committee 
of the Agricultural University of Athens under the guidelines of 
"Council Directive 86/609/EEC regarding the protection of ani- 
mals used for experimental and other scientific purposes". 

Adipose tissue protein extraction and protein content deter- 
mination prior to electrophoresis 

Adipose tissue frozen blocks were grinded in pestle and mor- 
tar with liquid nitrogen and the powder was collected and mixed 
with Radio-Immunoprecipitation Assay buffer (RIPA) buffer 
(R0278, Sigma-Aldrich) containing 1 mM Phenylmethanesulfo- 
nyl fluoride (PMSF) (P7626, Sigma-Aldrich) and protease inhib- 
itor cocktail (Roche Diagnostics, 11836170001). Lysates were 
concentrated by centrifugation and supernatants transferred in 
new pre-chilled tubes, snap frozen in liquid nitrogen and stored 
at -80 °C for further analysis. 

Total protein from the adipose tissue samples was determined 
by Bicinchoninic Acid (BCA) protein assay (23227, Thermo 
Fisher Scientific) in a 96-well plate according to manufacturer 
instructions. Bovine serum albumin (BSA) standard curve was 
constructed from various concentrations of BSA in RIPA buffer. 
The BCA protein assay reaction mixture added to the unknown 
protein sample as well as standards incubated at 37 °C for 30 min 
and immediately measured at 562 nm by using a microplate 
reader. The protein concentrations were determined by compari- 
son to the standard curve. 

Sodium dodecyl sulfate PAGE (SDS-PAGE) and immunob- 
lot analysis 

The expression of 0G6PDA and 0G6PDB was assessed in adi- 
pose tissue immunoblot analysis by using goat antibody raised 
against 0G6PDH. Prior to this investigation, in order to identify 
the presence of these variants, epitope peptide tracking was based 



on evidence found by in silico study, searching for the antibody 
immunogen sequence homologs in the protein database. The 
above study equally predicted the putative 0G6PDH protein iso- 
forms, and confirmed the capability of this antibody to capture 
equally the protein targets. The protein blast ([BLASTP 2.2.28] 
citation http://blast.ncbi.nlm.nih.gov/Blast.cgi) for the immu- 
nogen peptide used for the antibody (Acris Antibodies GmbH), 
identified among other hits for G6PDH, the variant A (acces- 
sion number: ABD34655.1), and variant B (accession number: 
NP_001087249.1) from Ovis aries. 

The technique was performed as previously described else- 
where. 35 Specifically, as a primary antibody goat polyclonal and 
affinity purified antibody against G6PDH (diluted to 1:12500, 
AP16595PU-N, Acris Antibodies GmbH), was used. As a sec- 
ondary antibody, a horseradish peroxidase (HRP) conjugated 
anti-goat antibody (diluted to 1:7500, ab674l, Abeam) was 
used. As a detection system, the enhanced chemiluminescence 
reagents (GE Healthcare) for immunoblotting, were employed 
as well. Consequently, PVDF membranes were re-incubated 
with mouse (3-actin antibody (diluted to 1:1000, MAB1501, 
CHEMICON International, Inc.), ensued by a horseradish 
peroxidase (HRP) conjugated anti-mouse antibody (diluted to 
1:10000, RPN2108, GE Healthcare). The transferred protein 
gel bands on polyvinylidene difluoride (PVDF) membrane, 
were visualized by chemiluminescent detection on a medi- 
cal X-ray film (Fuji safety film, RX 13 cm x 18 cm, RX1318, 
FUJI) and quantified by scanning densitometric analysis using 
the ImageJ analysis software. 36 Results are from two indepen- 
dent blots, with averaged values from the densitometry (« = 2), 
while the densitometric data were normalized to (3-actin signal 
intensity. 

0G6PDH enzyme assay 

One sample from each adipose tissue was frozen at -20 °C for 
enzyme assays for further analysis. For enzyme assays, about 500 
mg of adipose tissue was homogenized in 0.15 M KC1 at 4 °C 
and centrifuged at 4 °C and 6000 g for 45 min. Duplicate mea- 
surements of G6PDH (EC 1.1.1.49) activities were performed 
photometrically at 37 °C, according to the methods described 
previously. 10 

The estimated enzyme activity was expressed in nanokatal 
(nkat), which is the SI unit (transformation of 1 mol of substrate 
per second). 

Data analysis 

Data are presented as mean ± SEM. We used one way analysis 
of variance to compare the data, and with a significant variance, 
we performed Bonferroni post test. Analysis was also confirmed 
by unpaired t test or nonparametric test (Kruskal-Wallis) fol- 
lowed by the Dunn multiple comparison test. All hypothesis test- 
ing was two-sided and a P value of less than 0.05 (P < 0.05) was 
accepted as significant. Linear regression analysis for repeated 
measures was also used to determine protein expression changes 
with energy balance (MJ). Subsequently, a posthoc analysis was 
applied for intercept (b) comparisons in the equations of the 
first grade (Fig. 4). Moreover, in order to investigate further the 
energy balance effect in 0G6PDH protein levels, the protein 
expression data were divided into two groups according to ewes 
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energy balance: (1) group of ewes with negative energy balance 
(n = 14); and (2) group of ewes with positive energy balance (n = 
7). Analysis was performed using GraphPad Prism version 5-04 
for Windows (GraphPad Software). 
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